Topological Resonances in Scattering on Networks (Graphs) 
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We report on a hitherto unnoticed type of resonances occurring in scattering from networks 
(quantum graphs) which are due to the complex connectivity of the graph - its topology. We 
consider generic open graphs and show that any cycle leads to narrow resonances which do not 
fit in any of the prominent paradigms for narrow resonances (classical barriers, localization due to 
disorder, chaotic scattering). We call these resonances 'topological' to emphasize their origin in the 
non-trivial connectivity. Topological resonances have a clear and unique signature which is apparent 
in the statistics of the resonance parameters (such as e.g., the width, the delay time or the wave- 
function intensity in the graph). We discuss this phenomenon by providing analytical arguments 
supported by numerical simulation, and identify the features of the above distributions which depend 
on genuine topological quantities such as the length of the shortest cycle (girth). These signatures 
cannot be explained using any of the other paradigms for narrow resonances. Finally, we propose 
an experimental setting where the topological resonances could be demonstrated, and study the 
stability of the relevant distribution functions to moderate dissipation. 
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Narrow resonances are abundant in a large variety of 
physical systems, and their immense importance as in- 
dicators of long lived states led to an intensive study 
of their properties and origin. The weak effective cou- 
pling to the continuum which underlies their appearance 
can arise in various circumstances. Common mechanisms 
are the existence of potential or dynamical barriers such 
as weakly transmitting optical mirrors in optical Fabry- 
Perot resonators, dynamical tunneling in systems with a 
mixed phase space P, 0] or scarring of quantum states 
by unstable trapped orbits @. Anderson localization 
in disordered systems [1, [H induces narrow resonances 
not because of any barriers but because of destructive 
interference between multiply scattered waves with ran- 
dom phases while a classical particle diffuses unhindered 
through the system. Here, the resonance parameters de- 
pend on the particular realization of the disorder and 
must be studied with statistical methods. Statistical 
methods are also necessary for stud^ng resonances which 
characterize chaotic scattering |6l-[lll| , where fluctuations 
in the wave functions may lead to approximate bound 
states with very low amplitude at the interface between 
the interior of the system and the continuum. This in- 
hibits the transition to the exterior, resulting in a long 
living state. 

Wave propagation on bounded networks (graphs) dis- 
plays many features which arc typical to quantum chaotic 
systems [ij, [l^ ■ When the networks are connected to ex- 
ternal leads, the resulting scattering parameters fluctu- 
ate, much in the same way as expected from the analysis 
of chaotic wave-scattering in open Hamiltonian systems 
H, However, in addition to resonances from 
random-like interfering waves, the non-trivial connectiv- 
ity is also responsible for the formation of another type of 
narrow resonances - the subject of the present note. We 
will show that these resonances exist in a large class of 



scattering graphs if the graph contains a cycle. We thus 
call them 'topological resonances'. They have properties 
which clearly distinguish them from other mechanisms 
leading to narrow resonances: their mark on the distri- 
butions of the resonance parameters cannot be explained 
by any of the other paradigms. 

The rest of this letter is organized as follows. The topo- 
logical resonances signature will first be illustrated with 
some numerical simulations. The underlying theoretical 
framework, will then be reviewed, and used to derive the 
observed resonance distributions in simple cases. Finally, 
a possible experimental setup is proposed where topolog- 
ical resonances could be observed. 
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FIG. f . (Color online) Tails of the integrated distribution I{a) 
for the mean intensity a. Insets show the underlying graphs. 
Full lines represent numerical results where L = 1, . . . , 5 is the 
number of attached lead. Note that the data for L = 2, 3, 4, 5 
for the cube consists of more than one line - these correspond 
to non- equivalent ways to attach L leads. The dashed lines 
are power-law tails with indicated exponents. 

Numerical illustrations: The insets of Fig. [1] show two 
different networks: a fully connected graph with 4 ver- 
tices (tetrahedron) on the left and a graph with 8 vertices 
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(cube) on the right. The interior of the graph consists 
of finite bonds between vertex pairs (red hnes). Infinite 
leads (gray hncs) are attached to some vertices. Scalar 
waves propagate freely on the bonds and leads, and at the 
vertices they are reflected or transmitted without losses. 
Thus, one has a scattering system, which is described by 
'ATI L X L scattering matrix S{k). 

Solving the scattering problem numerically with an in- 
coming wave of unit flux and wave number fc, we com- 
pute the mean intensity a{k) of the wave function on the 
internal edges (see (jU for an explicit definition). The 
resulting distributions 



P{a) = lim — 



K 



5{a — a{k))dk 



(1) 



or rather, their cumulative form /(a) = P{a')da' , 
display extremal power-law distributions shown in in 
Fig. [TJ for various L values. The simulation summarized 
in Fig. [T] and further numerical results for other graphs 
suggest /(a) ^ for large a with 



for L < C - 1 



''=14^ forL>C-l 



(2) 



where C is the number of bonds on the shortest cycle 
(girth). The tetrahedron girth is C = 3. Thus, one 
expects the exponent to be = 1-5 for i = 1 and 
jjLL =2 for L > 2. On the other hand, for the cube with 
C = 4 one expects fii = 1.5, /i2 = 2, and fiL = 2.5 for 
L > 2. This is indeed borne out by the simulations. The 
appearance of the girth suggests a topological origin. 

The result ^ cannot be explained by the other 
paradigms for narrow resonances. Indeed barriers such 
as almost perfect mirrors in a Fabry-Perot interferome- 
ter lead to a cut-off for the intensity. Scattering from a 
disordered system gives rise to a power law with a fixed 
exponent uioc = 1 which is independent of the number of 
channels [J, Q . Random- matrix models for chaotic scat- 
tering predict a power law with an exponent ^rmt = 



which is only consistent with ^ if L < C — 1. 
We shall now summarize the graph theoretical setting 
(including the definition of the class of graphs) where 
topological resonances will be rigorously defined, and 
equation ^ will be derived. 

Quantum Graphs and topological resonances: A scat- 
tering graph G = {V,£} consists of a set of vertices V 
and a set of edges £ = BU C where B is the set of bonds 
connecting pairs of vertices and C the set of infinite leads. 
The graphs considered here arc all of finite cardinality. 
We will assume that the graph is connected, has no loops 
(each bond connects two different vertices) and each ver- 
tex i has degree (the number of attached edges) di > 2. 
Thus, graphs with dangling bonds corresponding to ver- 
tices with d = 1 are excluded. Multiple connections be- 
tween two vertices are allowed. Each bond b G B has 
a finite length £b € (0, oo) and a coordinate Xb £ [0, £b] 
such that Xb = and Xb = ib correspond to the two end 



vertices. The leads I € C are of infinite length. The co- 
ordinate xi £ [0, oo) is defined such that a:; — is the 
position of the end vertex. 

The complex valued wave function on the graph is writ- 
ten as * = {Mxb)}beB U {Mxi)}iGC = {V'e(2;e) legi- 
on each edge the wave function satisfies the Helmholtz 



(or stationary free Schrodinger) equation 



dx'i 



+ k^e = 



where fc > is the wave number. At the vertices the wave 
function is continuous and X)e=i ^^('^) ~ where the 
sum extends over all edges which emanate from the ver- 
tex. These Neumann or Kirchhoff matching conditions 
are a standard choice from a wider range of admissible 
boundary conditions [l6j . At a given wave number the 
wave function on any edge has the form 



-ikxe 



(3) 



where a(e,-) and a(e,-i-) are the amplitudes of the two 
counter-propagating waves on the edge. The mean in- 
tensity a{k) on the graph is defined by 



' ' beB;cr=± 



ab.o 



(4) 



If the bond lengths are rationally independent then 
the spectrum of the quantum graph is purely continuous 
with generalised eigenstates which are bounded every- 
where but are not normalizable. To each value of fc > 
one can associate a unitary scattering matrix S{k) 
that relates the outgoing coefficients on the leads to the 
incoming ones. Resonances are identified as poles of the 
scattering matrix when k is in the upper complex fc-plane. 
The (positive) imaginary part of the wave number at 
a resonance gives the decay rate (width). If the bond 
lengths are changed continuously so do the positions of 
the poles of S{k). When bond lengths become ratio- 
nally dependent some poles may move to the real axis 
indicating the appearance of a normalizable bound state 
embedded in the continuum (see [T7l - [l9| ). 

We can now define a topological resonance as a pole of 
the scattering matrix that can be moved to the real line 
to form a bound state by changing some bond lengths 
continuously while the graph connectivity and matching 
conditions remain unchanged. 

The main statement of this letter is that for rationally 
independent bond lengths a quantum graph as described 
above supports topological resonances if and only if it 
contains a cycle. 

This can be shown following similar ideas as in [20| . 
Let us first consider any scattering graph that does not 
contain a cycle. Such a graph is a tree (and to satisfy the 
requirement d > 2 all the canopy edges must be leads). 
The wave function for a bound state has to vanish on all 
the leads in order to be square-integrable. The remaining 
edges (bonds) form a finite tree whose leaves are at one 
end only connected to leads. The matching conditions 
then imply that the wave function also has to vanish on 
all the leaves. The process can now be iterated on the 
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remaining tree, showing that there are no non- vanishing 
square-integrable solutions on a scattering tree graph (of 
finite cardinality). As a result there cannot be any topo- 
logical resonances on a scattering graph without cycles. 
Now let us assume that the scattering graph contains a 
cycle and let us show that there are choices for the bond 
lengths that lead to bound states which we call topological 
bound states. Topological resonances are the remnants of 
topological bound states when the latter get mixed with 
the continuum by a generic choice of bond lengths. The 
mechanism is related to a similar phenomenon for closed 
graphs where it explains the structure of scars (20j . Let C 
be a cycle in the graph that consists of C = \C\ bonds as 
shown on the left in Fig. [21 Our assumptions (no loops) 
imply C > 2. Now let us construct a bound state on 
C. We require that the wave function vanishes exactly 
outside C. By continuity all vertices on C are then nodal 
points. The wave function for any b £ C then has to be of 
the form ±sin(A;a;ti) and the matching conditions reduce 
to the statement that the wave function and its deriva- 
tive have to be continuous along the cycle. This implies 
the following two diophantine conditions: 

kib X TT with nf, £ N for all 6 £ C (5) 
rib ~2s for some s G N. (6) 

bee 

The above conditions can be satisfied for a discrete se- 
quence of wave numbers if and only if the bond lengths 
on C are rationally dependent. This means that there 
exists a unit of length £o > such that £b = ib^o for some 
ib G N. One finds topological bound states on the cycle 
for wave numbers nkg (2n/co) with n G N if J2bee 
even (odd). The only reason for the existence of such 
states is the combination of a topological structure (the 
cycle) with destructive interference at the nodes. 
A generic (rationally independent) choice of bond lengths 
will destroy topological bound states on the cycle as the 
conditions (O and © cannot be satisfied exactly. How- 
ever the condition can be satisfied approximately to 
arbitrary precision at appropriate wave numbers. This 
is where the destroyed topological bound states leave a 
mark in the form of topological resonances. 

Let us now show that topological resonances on cy- 
cles indeed lead to the statistical signatures that we ob- 
served in the numerical simulations reported above. For 
this we consider a connected scattering graph with girth 
C > 2 and focus our attention on the corresponding cy- 
cle. Each vertex has two attached bonds which belong to 
the cycle and at least one additional attached edge that 
does not belong to the cycle (as shown in the left part of 
Fig. [5]). We will show that the ratio p of the intensity in- 
side the cycle to the intensity outside has a distribution 
with a power- law tail P{p) ^ p^^'^^^ with pc = 
For simplicity we will assume that all vertices on C have 
the same degree d ~ 3 (the calculation for the general 
case follows the same steps but is too cumbersome for 
this note). Since we only want to compare the intensity 
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FIG. 2. (Color online) Left: A cycle of length C. 
Right: Numerically obtained tail of the integrated distribu- 
tion I (a) of the mean intensity for a lossy beam-splitter setup 
(see inset). 

inside the cycle to the intensity on the edges that are ad- 
jacent to the vertices on the cycle we may disregard the 
rest of the graph and replace the adjacent edges by leads 
of infinite length. We have thus reduced the problem to 
finding the scattering solutions for one cycle of length C 
with one lead at each vertex. The mean intensity on the 
cycle is defined as in Q where the sum is restricted to 
bonds on C. Combining the amplitudes in one vector 
one finds (see [T^. ITsI) 

ac = (l - e*"/2acce*'='/') e'^'/^accK (7) 

Here € is a diagonal matrix of size 2C x 2C that con- 
tains each bond length of the cycle twice, b™ are am- 
plitudes of incoming waves on the leads such that the 
mean intensity outside the cycle is proportional to |b™p 
(flux conservation ensures that the outgoing waves have 
the same intensity). The 2C x 2C matrix acc and the 
2C X C matrix acc contain scattering amplitudes that 
can be derived from the Neumann matching conditions 
at the vertices. They are given by 

where P is the permutation matrix for the cyclic per- 
mutation (12... C). The sub-unitary matrix acc has 
one eigenvalue equal to one with the eigenvector vg = 

^ ^-|t ^ where Ic is the C-dimensional vector with unit 
entries. Bound states appear whenever 

^k) = e'*^-^/Vcce'"/' (9) 

has an eigenvalue equal to one. This happens if e**^^/^ = 1 
which cannot be satisfied for generic (rationally inde- 
pendent) bond lengths and A; > 0. However, in that 
case k e'*^^^^ is an ergodic fiow on a C-torus and 
gifc£/2 _ J defines a point on the torus that can be ap- 
proached to arbitrary precision as k increases [2]|. Defin- 
ing the ratio of intensities as p{k) = jacp/lb^p we 
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may derive a power law tail for P{p) = {6{p — p{k)))k 
by replacing the spectral average by a torus average 
P{p) = {2tt)^'~' J d'~''9S{a — a{9)) where we have replaced 
the 2C X 2C-niatrix e^''^ by a diagonal matrix e^^ that 
parameterizes the C-torus (each angle 9j appears twice). 
Focussing on the contribution from a small region around 
6 = where S(0) has an eigenvalue one which dominates 
the behaviour. Second-order perturbation then yields 



Pi0) 



/(2)(0) 



(10) 



where = ^ J2n=i ^" ^'^'^ On = On — 0. The functions 
f'^'^\Q) and g^^^O) are (explicitly known [g^) positive 
definite quadratic forms in the variables On ■ This implies 
P{p) - p-^'^-i with pc = ^- 

Coming back to the full graph that contains C as a sub- 
graph, note that the mean intensity a on the graph con- 
tains a contribution proportional to p from the cycle. 
As a consequence the tail of P{a) cannot decay faster 
than the tail of P{p) such that pc gives a lower bound 
for the exponent p. This is consistent with our numeri- 
cal findings (0) if i > C - 1. The exponent p = for 
L < C — 1 is consistent with the random- matrix approach 
and can be derived as a lower bound for the exponent in 
the present context following similar ideas [Q-Q- Note 
however that our simulations were performed on regular 
structures where the connectivity does not vary strongly 
and where at most one lead was connected to one vertex. 
If we want ([2]) to be true in other circumstances we need 
to redefine L appropriately. E.g. one may show that if 
there are many leads at the same vertex with Neumann 
matching conditions there is still only one open quantum 
channel that couples to the inside of the graph. Moreover 
if a graph has a (possibly large) subgraph that is weakly 
connected to the rest (e.g. via a single bridge) then there 
is only a small number of channels which connect to the 
subgraph. The following definition will take care of these 
issues. Consider a connected subgraph H. A vertex v is 
on the boundary dH of that subgraph if it is adjacent to 
at least one edge in H and at least one edge outside H. 
We redefine L as the minimum of the size (cardinality) of 
dH over all connected subgraphs H that contain a cycle 
and that contain no lead. This reduces to the number 
of attached leads for the numerical simulations presented 
before. With this definition we conjecture that the re- 
sult (m we have found remains true for generic scatter- 
ing graphs with Neumann matching conditions. We have 
excluded loops as they always lead to topological bound 
states that cannot be destroyed by changing the length 
of the loop. One may allow loops if one defines C as 
the shortest cycle which is not a loop. Dangling bonds 
have been excluded because they lead to a further set 
of topological resonances on paths between two vertices 



with degree one (i.e. between two vertices with mirror- 
like reflection) - ^ will remain correct if C is redefined 
appropriately [2^ . Let us also mention that measuring 
the intensity on one point of the graph will in general give 
a different but predictable power-law exponent - they are 
related to the local topological structure (e.g. the small- 
est cycle that contains the point) rather than the global 
topology. Only very regular graphs like the ones we used 
in Fig. [1] have the same power-law exponent at every 
point in the graph (and thus also for the mean inten- 
sity). However, other global quantities which are used to 
characterize resonances, such as e.g., the Wigner delay 
time or the resonance widths can be shown to distribute 
as the parameter a which was used here to render the 
theoretical discussion more transparent. 

The derivation of the exponents ^ made explicit use 
of some properties of the Neumann matching conditions. 
The derivation can be generalised to other continuous 
matching conditions. For non-continuous matching con- 
ditions topological resonances as defined above can be 
constructed but more complex topological features are re- 
flected in the exponents of the corresponding power laws. 
In Fig. [2] (right panel) we give numerical evidence of a 
topological resonance in a graph structure with match- 
ing conditions that can be realized in a laser experiment 
using beam-splitters and mirrors. The numerics show 
a clear power-law distribution with exponent p = 3/2 
(dashed line). We have also included loss (e.g. at the 
reflection from mirrors) - we characterize the loss by the 
parameter A which gives the fraction of photons which 
are lost when travelling once through the whole system. 
Lossy setups follow the power-law behaviour up to a cut- 
off that increases as losses decrease. Reducing losses to 
A = 10~^ may be within reach [2^ showing that topo- 
logical resonances may be observed in experiment. 

The narrow resonances in networks play a very signif- 
icant role in the presence of a nonlinearity as present in 
nonlinear optical wave guides or active optical fibres. The 
enhanced intensity at a topological resonance amplifies 
the nonlinearity to such an extent that the perturbative 
treatment breaks down and typical nonlinear effect such 
hysteresis appear [23|. A detailed discussion and classi- 
fication of topological resonances is now in preparation 

[H. 
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